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Abstract 



The hehcity amphtudes of the high-energy photon sphtting in the external Coulomb 
field are obtained exactly in the parameter Za. The cases of screened and unscreened 
potentials are investigated. The consideration is based on the quasiclassical approach, 
valid for small angles between all photon momenta. New representation of the qua- 
siclassical electron Green function is exploited. General expressions obtained are 
analyzed in detail for the case of large transverse momenta of both final photons 
compared to the electron mass. 
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1 Introduction 



It is well known, that the virtual electron-positron pair creation in an external Coulomb 
field gives rise to such nonlinear QED phenomena as Delbriick scattering ( coherent photon 
scattering [jl|) and the splitting of one photon into two. At present the process of Delbriick 
scattering has been studied in detail both theoretically and experimentally (see recent 
review 0). At high photon energies uj ^ m {m is the electron mass, h = c = 1) the 
scattering amplitude has been found exactly in the parameter Za ( Z\e\ is the nucleus 
charge, a = e^/4n = 1/137 is the fine-structure constant, e is the electron charge). The 
approaches used essentially depended on the momentum transfer A = |k2 — ki| (ki, k2 
being the momenta of the initial and final photons, respectively). The main contribution 
to the total cross section of Delbriick scattering comes from small momentum transfers 
A ^ u; (the scattering angle ^ ~ A/uj ^1). In this case the amplitudes have been 
found in p, ^ ^ by summing in a definite approximation the diagrams of perturbation 
theory with respect to the interaction with the Coulomb field, and also in ^ within the 
quasiclassical approach. It turned out that a.t uj ^ m and Za ~ 1 the exact in Za result 
significantly differs from that obtained in the lowest order of the perturbation theory. 

The applicability of the quasiclassical approximation is based on the fact that, accord- 
ing to the uncertainty relation, the characteristic impact parameter is p ~ 1/A and the 
corresponding angular momentum is/~co'p~c<j/A^lat small scattering angles. This 
fact was used in Refs. @, 0, where the quasiclassical Green function has been derived from 
the integral representation for the Green function of the Dirac equation in the Coulomb 
field P] . In Refs. [^, ^ the quasiclassical Green function of an electron for an arbitrary de- 
creasing spherically symmetric potential has been obtained, which allowed one to calculate 
the Delbriick scattering amplitudes in a screened Coulomb potential. 

So far the process of photon splitting has not been observed, although some events in 
the experiment performed at DESY [|ri|] were erroneously interpreted as photon splitting. 
As it was shown in [|T^ these events were due to the electron-positron pair production 
accompanied by hard-photon bremsstrahlung. Some possibilities to observe photon split- 
ting have been discussed in |jl3|. Photon splitting has been investigated theoretically in 



1^, O, |T6|, |T3] in the lowest order of perturbation theory with respect to the parameter 



Za. The expressions obtained in [0, |15| are rather cumbersome and it is difficult to use 
them for numerical calculations. Nevertheless, some calculations based on the results of 
Ref. |]14|, ITH] have been carried out in ||l^ |18| . Using the Weizsacker- Williams method 



providing the logarithmic accuracy the amplitudes of the process have been derived in an 



essentially simpler form in Ref. |16]. The comparison of the exact cross section [17| with 
the approximate result |TB| has shown that at high photon energy the accuracy is better 
than 20%. The magnitude of the Coulomb corrections to the lowest order amplitude of 
photon splitting has been unknown up to now. At the present time the experiment dedi- 
cated to the observation of high-energy photon splitting {u ^ m) in a strong Coulomb field 
is held in the Budker Institute of Nuclear Physics (Novosibirsk). Therefore, the theoretical 
investigation of the problem is of great interest. 

In the present paper the high-energy photon-splitting amplitude is calculated exactly 
in Za for small angles /2 and /s between the momenta k2 , ks of the final photons and 
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the momentum ki of the initial one. It is the region of small angles that gives the main 
contribution to the total cross section of the process. Besides, small angles and high energies 
of photons allow one to use the quasiclassical approach developed in |, [1^ at the 

consideration of Delbriick scattering. We consider the case of a pure Coulomb potential 
as well as the influence of screening. The initial representation for the splitting amplitude 
is rather cumbersome and contains a thirteen-fold integral. The quasiclassical approach 
gives the transparent picture of the phenomenon and allows one to determine the region 
of integration which gives the main contribution to the amplitude. Without that it seems 
impossible to calculate the amplitude. 

Our paper is organized as follows. In Sec. II we perform some transformations of the ex- 
act amplitude. These transformations essentially simplify further calculations. Section III 
contains the discussion of the process kinematics. In Sec. IV the small-angle approximation 
for the quasiclassical Green function is derived. In Sec. V this Green function is applied to 
the calculation of the photon-splitting amplitude. In Sec. VI we consider the case of large 
transverse momenta of both final photons (^2/2 ^ 1^3/3 ^ The limiting cases of 
small momentum transfers and Za <^ 1 for this amplitude are studied in Sec. VII and Sec. 
VIII, respectively. In the last Section the cross sections obtained in the Born approximation 
are presented and compared to those obtained within the Weizsacker- Williams method. 



2 Transformation of the amplitude 

According to the Feynman rules the photon-splitting amplitude in the Furry representation 
can be written as: 

M = ie^ J d^xTr < xleie-'^i^ge^e'^^^^^^e^e'^^'^lx > +{k>^ ^ ^ e^} . (1) 

Here and 62 3 are the polarization vectors of the initial and final photons, respectively, e = 
e^7^ = —67, 7^ being the Dirac matrices, Q = 1/(V — m + iO), and = idfj_ + g^Q^Za/r). 
The matrix element of the operator Q is the Green function of the Dirac equation in the 
Coulomb field: G{x,x') =< x\Q\x' >. 

It is convenient to rewrite the expression (0) in the form, containing only the Green 
functions of the 'squared' Dirac equation: 

D{x,x') =< x\V\x >=< x\l/{V^ -m^ + iO)\x' > 

For this purpose we represent the left operator Q in (|l|) in the form Q = V{V + m) and 
use the commutative relation 

{V + m)e e'^'' = e'^"" [ -e Q'^ + eA; - 2ep ]. 

One can obtain another expression by similar transformation of the right operator ^ in (|lD. 
Taking a half-sum of these two expressions and using the identity ^ dxTi < x\AiA2\x > = 
J dxTr < x\A2Ai\x >, valid for arbitrary operators Ai and A2, we get: 
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e^Q*^Qi{k2-k,)x ^ x\Ve''''%elh - 2e*p)P|x > + (2) 



Q*^Q*^Q-iik2+ki)x ^ x\Ve-'^^''{-eiki - 2eip)P|s > 



<x|e-^'=^^(-eiA;i - 2eip)Ve'^^'' {e^h - 2e;p)Ve''''%elk3 - 2e;p)V\x> +(2 ^ 3) 



+ 



The formula represents the amphtude as a sum of the terms containing either two 
or three Green functions: M = M^^^ + M^^\ The amphtude in the form of eq. (|I]) contains 
terms of different order of magnitude and a strong compensation occurs. After taking the 
trace in eq. (0) the expression obtained contains only the terms of the required order what 
is the advantage of this formula. 

Passing from the time- dependent Green functions to the energy-dependent ones, taking 
the integral over time in eq @) and omitting the conventional factor 2Tr6{uJi — uj^ — ciJa), we 
get for the contribution to M containing three Green functions 

M(3) = -e^ y — y rfrirfr2rfr3exp[i(kiri - ksrs - kara)] x (3) 
Trj [i-eih - 2eip)D(ri, \€ - cua)] [{elh - 2e;p)D(r2, \e)] x 
[{elk - 2e*p)D(r3, ri \e + u-s)]] + {k^ ^ k^ , 62 ^ 63) . 



Here p = — iV differentiates the Green function D of the 'squared' Dirac equation with 
respect to its first argument. 

Below the Green function Z}(ri,r2 {e) will be called the 'electron Green function' for 
e > and the 'positron Green function' for e < 0. Let the initial photon propagates along 
z axis. Then according to the quasiclassical approach developed in p, ^ [1^ the main 
contribution to the amplitude at high energy arises from the region of integration over 
the variables Zi in which z' < z for the electron Green function D{r,r'\e) {e > ) and 
z' > z for the positron Green function {e < ). In terms of noncovariant perturbation 
theory this range of variables corresponds to the contribution of intermediate states for 
which the difference between their energy En and the energy of the initial state Eq = ui 
is small compared to Eq. Outside this region at least for one of the intermediate states 
\En — Eq\ ^ Eq and therefore the corresponding contribution is suppressed. Besides, there 
exist another restriction for the region giving the main contribution to the amplitude. This 
restriction follows from the explicit form of the quasiclassical Green function. It reads 

zi < Z2,Z3 , zi <0 , max(2;2, 2:3) > . 

All these conditions allow one to depict the main contribution to the amplitude M*^^-* in 
the form of diagrams, shown in Fig. 1. The explicit form of vertices is obvious from eq. 
(H). The electron Green functions are marked with left-to-right arrows, and positron ones 
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with right-to-left arrows. The arrangement of the diagram vertices is space ordered. With 
the use of these diagrams one can easily determine the limits of integration over the energy 
and coordinates. The diagrams (a) and (b) correspond to the following picture: the photon 
with momentum ki produces at the point ri a pair of virtual particles which is transformed 
in the point r2 into a photon with momentum k2. Between these two events the electron 
(a) or the positron (b) emits a photon with the momentum ks at the point r^. 



Analogously, the expression for the term M*^^) containing two Green functions reads 

de 



I ^ I ^^^^^2 Tr|exp[i(kiri - k2r2 - k3r2)] e^e^ x 



(4) 



[(-eifci - 2eip)D(ri, T2 \e - ui)] D{r2, ri \e)] 



exp[z(kiri - k2r2 - kgri)] eie* x 



D{ri, ra \e - UJ2)] [{elh - 2e;p)D(r2, ri \e)] + {k^ , e2 63) 



The diagrams, corresponding to the representation (H), are shown in Fig. 2. 



3 Kinematics of the process 

According to the uncertainty relation the lifetime of the virtual electron-positron pair is 
r ~ |r2 — ri| ~ ojilim? + A^), where A = max(|k2±|, |k3j_|) ^ uji , k2± and ka^ being the 
transverse components of the final photon momenta. The characteristic transverse distance 
between the virtual particles can be estimated as (m^-|-A^)~^/^, which is much smaller than 
the length of the electron-positron loop. The characteristic impact parameter is p ~ 1/A , 
where A = k2 -|- ka — ki is the momentum transfer. At small k2_L and ka^ (/2,3 ^ 1) we 
have 

A^= (k2, + ka,)^ + ^f^ + ^V. (5) 

A\ UJ2 UJ-i j 

The characteristic angular momentum is / ~ cj/A ^ 1 , and the quasiclassical approxima- 
tion can be applied. 

Let us discuss a screened Coulomb potential. In the Thomas-Fermi model the screening 
radius is Tc ~ {jnQ)~^Z~^l'^ . If i? <^ 1/A ^ Tc (-R is the nucleus radius), then the 
screening is inessential and the amplitude coincides with that in the pure Coulomb field. 
At 1/A ~ Tc the screening should be taken into account. Obviously, the impact parameters 
p ^ Tc do not contribute to the total cross section. Due to this fact we shall concentrate 
ourselves on the momentum transfer region corresponding to the impact parameter p < Tc- 
If k2_L ^/co'2 -|- ka_L ^/t^3 ^ T'^^ , then it follows from that the condition p <rc holds only 
when I A_|_| = |k2_L + ka^l > r~^. Thus, the main contribution to the amplitude is given by 
the region of momentum transfer A_|_, restricted from below. In this region |A_|_| ^ I'^yl) 
that is A ^ Aj^. It is this region of parameters which we are going to consider. In addition, 
at uj/{m? + A^) ^ Tf. the angles between ki 2,a and ri 2,a are either small or close to vr , and 
corresponding expansions are used in our calculations. 

According to the Furry theorem the photon-splitting amplitude is an odd function with 
respect to the parameter Za. Due to the singularity of the Coulomb potential in the 
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momentum space (— 47rZa/A^) the region of very small momentum transfers A < 
is essential only in the lowest (linear in Za) order of the perturbation theory . In this 
order the Weizsacker- Williams method is applicable and the corresponding large logarithm 
appears in the cross section integrated over the angles of one of the final photons ||T6[. In 
next orders of perturbation theory with respect to the parameter Za the integral over all 
momenta corresponding to the external field should be taken provided that their sum is 
equal to A. Therefore, even at A ~ each momentum is not small and the screening 
can be neglected. In the Born approximation the screening can be taken into account by 
multiplying the amplitude by the factor [1 — F(A^)], where F(A^) is the atomic electron 
form factor. Thus, to find the photon-splitting amplitude in a screened Coulomb field it is 
sufficient to solve the problem in the pure Coulomb field. 



4 Green function 

Let us pass now to the consideration of the Green function D{r,r'\e) appearing in eqs. 



(0) and (^. In Refs. 10] the representation of this function has been found in the 
quasiclassical approximation for an arbitrary decreasing spherically symmetric potential. 
For the case of the Coulomb field, at small angle 6 between vectors r and — r' we obtain 
from eq. (14) of g: 



D(r,r'\e] 



dll 



W 



X 



(6) 



exp 



J2(r + r') 
2nrr' 



Ak rr' 



where a = 7^7 , = 5^ — , X = e/n , n = r/r and n' = r'/r'. Taking into account 
the relations 



/ 



dllUie)g{f) = rfqexp(^q0)^(g2) , ^,(16) = -1 A 



where gil"^) is an arbitrary function and q is a two-dimensional vector, one can rewrite eq. 
(E) in the form 



D(r,r'\e) 



iK{r+r') 



1 dq 


1 + Za^r- 







(7) 



exp 



Iurr 



TV 



Here 6 = r±/r , 0' = r'j_/r'. Expression (0) contains only elementary functions and the 
angles and 6 ' appear only in the factor exp[zq (0 + ')]. Therefore, the representation (|^) 
for the Green function is very convenient for calculations. At small angle between vectors 
r and r', for the case of pure Coulomb field, we get from eq. (15) of H: 



D{v,Y'\e) 



47r|r 



^ ^ »ZQ:ASlgn(r-r') 



(8) 
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One can see that in this case the Green function differs from that of free Dirac equation 
only by the phase factor. It is easy to check that after the substitution of the expressions 

(0) and (H) for the Green functions to the sphtting amphtudes @ and (H) all phase factors 
of the form r^*^" cancel. Note that one can replace k by \e\ — m^/(2|e|) in eqs. (|^) and 

(1) . Moreover, one should take the relativistic correction m^/(2|e|) into account only in 
the factor exp[m(r + r')]. 



5 Calculation of the amplitudes M^^) and M^^) 

Consider now the diagrams, containing three Green function (See Fig. 1). Obviously, 
the contribution of diagram (b) can be obtained from the contribution of diagram (a) by 
replacing Za —Za and changing the overall sign. It provides the fulfillment of the 
Furry theorem: the sum of two contributions (a) and (b) in Fig. 1 is odd function of 
Za. Therefore, this sum is equal to the odd in Za part of diagram (a) multiplied by two. 
Then, at the calculation of diagram (a) the region of integration over is divided into 
two: Z3 > (photon with the momentum ka is ahead) and 2:3 < (photon is behind). In 
the first region the angles between vectors r2 , and — ri are small. In the second region 
the angles between ri, and — r2 are small. Denote the contribution of the first region 
to the diagram (a) as M^-* and the contribution of the second region as M2^'^ . Introduce 
vectors f2 = k2_L/t^2 and fs = k3_L/co'3 (|f2,3| < 1), as well as di = Yi^/vi = ni± (i=l,2,3). 
Taking into account the smallness of the angles we have dvi = rf dri dOi . It is convenient 
to perform further calculations in terms of the helicity amplitudes MAiA2A3(ki, k2, k3). It 

is sufficient to calculate three of them, for instance, M_| (ki,k2,k3) , M+++(ki, k2, k3) 

and M++_(ki, k2, k3) . The others can be obtained by substitutions. Note that within 
the accuracy of our calculations there is no need to take into account the corrections 
to the transverse part of the polarization vectors 62,3. Due to the relation ek = the 
z— component of the polarization vectors can be expressed via the transverse one as (62,3)2 = 
— e2,3f2,3 • So, for a given helicity one can put the transverse part of the final photon 
polarization vector to be equal to the polarization vector of the photon with the same 
helicity propagating along z axis. Further the polarization vector e+ corresponding to the 
positive helicity is denoted as e. Then the negative helicity polarization vector e_ is equal 
to e*. Note that to obtain M++_ in our approach it is necessary to calculate both M\i__ 
and MXL+ since the arrangement of the vertices on the diagrams is space ordered. 

Let us substitute the expressions (|^) and (|^) to (^) and make the obvious expansion 
at small angles, taking into account terms, quadratic in fj and 6i. Introduce notations 
Hi2 = UJ2 — e , K3 = u;3 + £ and pass to the variables 

D ^1 D ^2 D ^3 

K2H3 SH2 SK3 



Simple integration over 6i leads to 

3 1^2 ooooL , . 2iZa 
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where L = i?2^^3't2/i^2%; and 



T 



Tr 1 



Ri 



ei Q - eih 1 



R 



R 



'el,Q + eR2U3)-elk: 



Za Q:q2\ 

Ai 



(10) 



R 



£^R2Rz f 23 

2 2i? 



[/t2q2 - Ksqa, A) 



(o^3^2i^2 - ^^2^:3^3) ir^f . I 

^ (Qf23) - — + R) 



R = Ri 



R3 



[23 



f2-f3 



Q = q2 + q3, A = t^2f2 + t^3f3 



It is convenient to transform the function T in eq. (|TOD to the form which does not contain 
the parameter Za. To this purpose we use the identities 



7 q2 / ?2 
Za^ — 



2iZa 



i d 



2 dq^ 



02 

.(13, 



2iZa 



13 



2iZa 



i d 



2 dq. 



02 
.(13, 



2iZa 



and integrate by parts over q2 and qg in After that some terms in the function T are 
independent of the variable Ri and others contain it in the form of factors 1/Ri, ^/R\- 
Taking the trace and integrating by parts over Ri the terms with the factor I/-R1, we obtain 
for different polarizations: 



[eQ)(eQ2)(eQ, 



3; ) 



1^2 ^ 1^3 



RiR^ Vk3 K2 



(eQ)(e*Q2)(e*Q3) - 



uj2_ 
sR^ 



e , Q3Q2 Q2Q3 



1^2 



K-3 



+ 



n2K2,R^ 

4 



[e*, Q2 + Q3) + 



eR 



RiR^ \e K2 



e , 10^2 ^3 

1^2 1^3 



UJ2UJ3 



(eQ)(eQ2)(e*Q3) + —^(Q^ + Q2)(eQ) 



to'iu;2 



Q 



K2K3R \ 



^-m'\ (eQ2) + 



4i / oJiOJi 



R\R \ Kj2^3 



-2 (eQ) 



2iujiUJ2 
K2K3-R2 



e, Q2 + Q 



3; ) 



where Q2 = Q + £-R2f23 and Q3 = Q + e-R3f23- The function T+_+ can be obtained from 
T++_ by the substitution c<j2 ^ 0^3 , K2 K3 , Q2 ^ Q3 and e —e. 



(3) 

In the same way for the contribution we obtain: 

„ UJ2 00 00 Li 



(3) 



?>2'k'^UJiUJ2^^3 



£K2K^de / dRi / (ii?2 









where Li = RiU3K2/uJie , r = Ri + R3 , 







'dR-i 



dci2 dq^ 



02 

.(13, 



2iZa 



e**T. 



(12) 



1 1 \ kIR,R^{ 
r R2) 2 2r 

{uJzl^2Rl - SUJiR;^) 



[/t2q2 -gq3, A) 
UJ2 



+ 



0J2r 



(Qf3)-— (/?2 + r) 



(13) 
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and the function T for different polarizations reads 



r2R 



■(eQ)(ePi)(eP3); 



(14) 



= ^ (- + - ) (e*Q)(eP0(eP3; 



^e,^P^P, + ^^P?P3 



K2er 



2 V"'' 



Kar V K2 



' ^^ + ^] (e*Q)(e*P,)(eP3) + -^(P? + PD(e*Q) - 

2 (eQ)-^— ^(e,Pi + P3; 



r2i?2 \K3 K2J 



K2er \ / 



rR2 \ K2e 



K2er^ 



where Pi = Q + K3i?if3 and P3 = Q — K3i?3f3. One can get the function from 
by the substitution ui ^ —u^ , /t2 ^ — e , Pi P3 and e ^ e*. Note that the integrand 
in eqs. (0) and ([I^ for hehcity amphtudes M2^^ can be obtained from the integrand in 
eqs. (D , (InD for Mf ^ by means of substitutions 



^2,3 -^2,3 > UJl^ iU2 , C^3 ^ "^^3 , /«3 ^ ^ 
i?l <-> i?2 , -R3 ^ — -R3 5 f23 ^ — f3 5 f 2 — f 2 



(15) 



so, that 

T+__^T+— , T+_+^T+-+, r+++ ^ r++_(e <-> e*) , r++_ T+++(e ^ e*) . 
As to the amphtude M^'^\ after the integration over the angles 6i, we get: 



m|'2_ = , MS+ = (e%Mi2 + Mi3) 



(2) 



(e,Mi2 + M23) 



(2) 



(e , Mi3 + M23) 



(16) 



where 
M23 



Ll)2 00 00 



de 



dRi f dR' 



Wn-' J ■ J RfJ Ri 

-LOS ^0 ^ 



Ri + 



R2 



dq^ dq^Qi — 



2iZa 



X 



exp< t 



1 1 \ Q\ C^2t^3«:2/t3n2 D ('«2q2 - '«3q3,^) ^ 



i?l ^ i?2/ 2 ^ 



2uj 



2 f 23-^2 

1 



-—{R, + R,) 



1^2 00 00 



26" /■ /-c/i^i /■rfi?2 

Mi3 = T7-^ 

^ 



167r3 



002 



d^2 dqsQ 



02 



2iZa 



exp< t 



1 , 1 \ U;it^3£:K2 2d ('«2q2 - £^q3> ^) 

^^^^^^ " - 



Ri R2 ) 2 



a;2 



^(^1 + ^2) 



(17) 



and the vector M12 can be obtained from M13 by the substitutions UJ2 ^ uj^ and f3 ^ f2. 
As we shall see, many terms cancel out in the sum M^^-* + M^^^ . 



In general case further transformations of the formulae obtained leads to four-fold in- 
tegral with the integrand containing the elementary functions and this problem will be 
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considered in detail elsewhere. In what follows we restrict ourselves to the case |k2_L| = 
\uj2^2\ ^ TT^i |k3±| = l^sfsl S> m when the amplitudes can be essentially simplified. This 
range of the parameters corresponds to large virtuality of the electron-positron pair com- 
pared to the electron mass, which in this case can be neglected. Note that the ratio 
between the momentum transfer A = |A| and the electron mass m can be arbitrary, since 
A determines a typical impact parameter p ~ 1/ A rather than pair virtuality. 



6 Zero mass limit 

It is easy to see that putting m = in the expressions obtained leads to the logarithmic 
divergences in some terms (in other words, these terms contain logm at finite mass). For 
instance, one obtains such logarithm in the amplitude M{ integrating over Ri the terms 
in T which do not contain the factor 1/Ri (see ([Tl|) ). The final result, as it should be, 
does not contain the logarithms of mass. But the cancellation of these logarithms between 
different terms is rather tricky. 

Taking the integral over Ri in (^ for T = we do not obtain any logarithm in 

M^\_ while M+l_ = 0. Eq. (|ll]) for T+_|_+ and T++_ contains the terms proportional 
to Q2. It is convenient to pass in these terms from the variables R2 and R^ to R2 and 
y = R3/R, and then integrate by parts over y. In the terms containing Q3 we pass to the 
variables R3 and y = R3/R and also integrate by parts over y. As a result, in the two-fold 
integral over R2 and R3 all logarithms cancel out and one can put m = 0. After that T^++ 
and T++_ are transformed to 

r+++ = -;^(^ + ^)(eQ)(e-Q,)(e-Q3) ; (18) 

In addition, there are integrated terms at y = u!3K2/^i£ ( upper limit) and y = (lower 
limit ). Making the similar transformations for the amplitude M2^^ we obtain that 
and T_|_+_ turn to 

T++- = ^(^ + -)(e*Q)(ePi)(eP3) ; (19) 

T 1X2 \ ^ ^2 ' 

The integrated terms corresponding to the lower hmit cancel out in the sum of Mf^ and 

M^^ . Remind that to calculate the amplitude M*^^^ we have to find the sum Mf ^ + (k2 ^ 

ks , 62 *^ 63) , extract the odd in Za part and multiply the result by two. After that 

(3) 

the contribution of the upper-limit integrated terms of M\ vanishes in the case 
and in the case M++_ and M+_+ gives the finite result after the summation with (eM23). 
To cancel the logarithmic terms we exploited the antisymmetry of some integrands with 
respect to the substitution e 002 — 'jJj, — ^ , q2 ^ "Qs- The contribution of the upper-limit 
integrated terms in M2 vanishes in the case of M++_, and for and M+_+ gives 
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the finite result at m = after tlie summation witli (e*Mi3) and (eMia), respectively. 
Analogously, the amplitudes (e*Mi2) and (eMi2) cancel singular terms of M2 (k2 ^ ka ) 
for the amplitudes and M++_. For the amplitude M+_+ the upper-limit integrated 

terms of M2^^(k2 ^ ks ) cancel out. As a result, the sum of integrated terms and M*^^-* 
gives the additional contributions to helicity amplitudes. We represent them in the form: 



5M 



47^3 



dR 



2iZa 



2iZa 



(20) 



For different polarizations the function F reads 



, = (eQ) 


I 






/ de „ ^ 









(ie — ^e" 



, f 2 ^ 1 

= -F+-+(u;2 ^ , f2 ^ fs) ; 

(K2q2 - KsQs, A) 



(21) 



Q CJ2U;3«:2/t3„2 



2a;2 



23 



The phase ip2 can be obtained from ipi by substitution (|15|). To obtain eq. ( |20| ) we have 
integrated over one of the radii. Note that the singularity of the integrand in ( pT]) at e = 
disappears in the total expression for the amplitude of the process so that no regularization 
is required. 

Putting m = in phases $ and $, we take elementary integrals over Ri in (|^) and over 
i?2 in (HID- After passing from the variables q2 and ga to Q = q2 + and q = q2 — q3, the 
integral with respect to q reads: 



J 




2iZa 



exp( 



:qA) 



(22) 



Let us change the variable q — |Q|q and represent Q and A as Q = |Q| Ai and A = | A| A2, 
respectively. After that it is easy to see that J depends on 5* = (A1A2) and |Q||A|. Note 
that in two-dimensional case the form P = eijX\X2 is also invariant with respect to rotations. 
However, the even powers of P can be expressed via 5* (P^ = 1 — 5*^), and the odd powers 
of P change their sign after reflection. On the other hand, J is invariant under reflection. 
It becomes obvious if the reflection of q is made along with the reflection of Ai and A2. 
Therefore, J is invariant with respect to the change of variables Q, ^ A.. Therefore, we 
can represent J as 

rfq /|q+ A| 
A^ ^"P^ 

This form is very convenient for further calculations. Now one can easily take the integrals 
over Q and all radii. Summing all contributions we finally get: 



2iZa 



J= ^ 



2^Q) 



(23) 



M 



8e^ 



7r'^UiU2UJ3A.' 



J dq(TVq)Im 1^ 



2iZa 



(24) 
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tJ2 

^ 



(eb)K3 (ec)e 



+ 



'UJ2 ^ 

V f 2 ^ f 3 y 



UJ3 J de K2 




UJ2'D-i 

e*b)(/ti + /t^) 
2(ea)cJiPi 



e f ti3-t^2 , (e*f3), 2 , 2 A 

^2 + '^sJ 



+ 



^UJ2 ^ UJ3\ 



LU2 

UJ3 de K,2 







1^3 



uJiVi 







a;2 e / K3 



(e*b) 



' 1^2 - £ _ (ef23) 
(ea) 

'e*f2 



2 

(e*f23) 



(ec)(/t^ + £^) 
2(e*a)cu2r'3 



+ 



£^3 + kI e"^ - ^2^3 

2^^11)1 2c<;3r'2 



+ 



Obviously, T+__|_ can be obtained from T^^.. by the substitutions UJ2 ^ uj^, ¥2 ^ fs- In 



eq. 



the following notation is introduced: 



UJ2K3a^ - t^3/t2b2 _ UJiKsC"^ - UJiSh^ 

Vi = tO , V2 = , V3 



UJlE UJ3K2 

a = q - A + 2K2f 2 , b = q + A - 2K3f 3 , 
c = q+A-2rf23 , c = q-A + 2£f23 • 



UJ2K3 



(25) 



At the derivation of (^41) we used the identity 

Qexp(--qQ) = 2«Vqexp(--qQ) 

and integrated by parts over q. Note that vectors e and e* appeared in denominators in 
(0) owing to the application of the relation 2(ea)(e*a) = a^. 



7 Asymptotics at small A . 

In the small-angle approximation (|f2|, |f3| ^ 1) the cross section of the process reads: 

^2 

da= —^\M\^x{l-x)dxdUdi3 , (26) 

where x = UJ2/UJ1, so that = uJi{l — x). Let us define p = {^2^2 — ^^3^3)/'^- In terms of 
the variables p and A the cross section has the form 

,,^,n dAdpdx 

Consider the asymtotics of the amplitudes at |A| ^ |p|. To this purpose multiply T in 
(0)by 

i = ^(go-q')+^(q'-9o) , 
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where |A| -C go |p|- Then, for the term in (|2^) proportional to ^{qQ — q^) one can 
put q = and A = in T and integrate by parts over q. After that, using the relation 
Vq'(9(gQ— q^) = — 2q(5(gQ— q^) one can easily take the integral over q since at |q| = go ^ 1^1 
one has 



Im 



|q+ A| 
|q - A| 



2iZa 



4:Za 



qA 

q2 



As a result, in the region |q| < go the term proportional to Za is independent of go and 
the terms of next orders in Za are small in the parameter |A|/go. 



For the term proportional to {}{a^ — q^) we get 

2iZa 2 



Vqlm 



AZa 



|q- 



q^A - 2q(qA) 



We put A = in T and perform the integration first over the angles of q and then over its 
modulus. As a result, the main in go/|p| contribution is independent of go and proportional 
to Za. Taking the sum of the contributions from these two regions and performing the 
integration over the energy e, we get 



4:iN{epf 



[Axp) 



N 



AZae^uj2UJ3 



N 



e*A + 2(eA) 



l + "^-^^ ln^ + 



UJi 



U2 



^2 + ^3 '^3 
2u'l UJ2 



TT 



N 



eA + 2(e*A) 



1 H fin h i-K) 

LU2 



^1 + ^3 n 2 ^3 



-(In^ 



+ 2m In 



t^3> 



2UJ2 ' OJi UJi 



where is the projection of the vector A on the direction of ki. Substituting (|28|) into 
(p7|) and performing the elementary integration over the angles of vectors A and p, we 
come to the expression 

AZ'^a^ dp^ d^^ dx , , 



da 



where the function g{x) for different polarizations has the form 

5f+__(x) = x(l - x) , 



(29) 



(30) 



-x(l — x) 



1 + 1 + (2a; - 1) In 



X 



X 



+ 



x^ + (1 — x) 



In' 



1 — X 



X 



1 , 
-x(l 

2 ^ 



x] 



1 + (2/x- l)(ln(l -x) +Z7r) + 



1 + (1 



x] 



2x2 



-(ln2(l -x) + 2i7r In(l-x)) 



9+A 



X 



Formulae (|29| ) and (^) are in agreement with the corresponding results of [|16] , obtained in 
the Weizsacker- Williams approximation. However, this approach does not allow to obtain 
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the amplitudes (28) themselves. The large logarithm appears after the integration of (^) 
over from A^j„ up to where A^m ~ '^c^ for the screened Coulomb potential and 
^min ~ 1^1 for the pure Coulomb case. It is interesting to compare the contributions of 
different helicity amplitudes to the cross section at A — » 0. In Fig. 3 the function g{x) is 
shown for different helicities as well as the quantity 



9{x) = g+—ix) + g+++ix) + ^++-(x) + g+-+ix) 



(31) 



which corresponds to the summation over the final photon polarizations. It is seen that 
g{x) has a wide plateau. 

The Coulomb corrections to the photon-splitting amplitude at A — > are small com- 
pared to the Born term (pHj). One can show that in this case the Coulomb corrections are 
proportional to (Za)^Aln^(A/p). Thus, they become essential at A ~ p. We shall present 
the detailed investigation of the Coulomb corrections elsewhere. 



8 Born approximation 



As it was mentioned above, the photon-splitting amplitude obtained in the lowest Born 



approximation [0, |T5| at arbitrary energies and momentum transfers is cumbersome and 
rather difficult to use. Therefore, it is interesting to consider the first term of Za expansion 
of (p4D. To obtain this asymptotics we perform the substitution 



(eVq) Im 



|q+ A| 
|q - A| 



2iZa 



Za 



(e*,q 



(e*,q 



in (|2^) and write the quantities Pi_3 in (|25|) as 



q+— A 

UJ2 



q + 



1^2 - 1^3 



f2 
^23 



iO 



+ 



f2 

^^3 



q 



/«3_+e 



(32) 



^Ux^2<^?,<^ f2 



After that we shift the variable of integration q in each term so that the quantities Pi_3 
become independent of the angle of the vector q. For instance, in the terms containing 
Pi we make a substitution q ^ q — As a result, one can pass to the variable 

z = exp(z0) and take easily the contour integral. Taking the integrals with respect to |q| 
and e, we get for the Born amplitudes 



2iZae\h X fg). 
7rA2(e*f2)(e*f3)(e*f23) 



(33) 



2{Za)e^uJi 



TtA (ef23)^Cc;2t^3 

Vr^ 1 , 2/'^2 



(eA) 



(et23) as 



(ef2)2 + (ef3)2 fn' , 1 



(ef 



23) 



, ln^(— ) + Li2(l-a2) + Li2(l-a3) 
6 2 as 



+ 
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feAl 



-^14 
1- a2 \ 

2(ef2)(ef3) 



02 ln(a2)\ 



(ef: 



23 J 



1-02 / 

02 In (02) 



+ CJ2(ef2 



a2 



002- 



Q3 

1 - as 

03 11 1(03) 

as 



1 + 



as ln(as) \ 
1-as / 



+ 



M. 



++- 



2{Za)e^uj2 
7rA^(e*fs)^co'icj3 



(e*A) 



(e*f2)2 + (e*f2s)V![! , lin^fZ^' 
V 6 2 ^ a2 ■ 



(e*A] 



^3(e*f23 



e*fs)2 
2 '^2 



(e%) + (e-f2s) , 
(e*f3) ^ a2 ^ ^ 

+ Li2(l -a2) + Li2(l + ai) I + 



1-^2 



'l + ^^MM^-c.?(e*f2)2 



ai 



1 - ct2 / 



1 + ai 



1 



2(e*f2)(e*f2s) 
(e*fs) 



ai ln(— ai) a2 ln(a2) 
iOi — — UJ^- 



1 + ai 



1 - a2 



ai ln(— ai) 

1 + ai 



+ 



where 



dt 



^ ^ ^ / N /■ 1 / N 

«1 = 72~ ' ^2 = 72 ' «3 = 72 ' Li2(x) = - / — ln(l -t) 

a;2ti^3i23 Ci;iu;2i2 tuitusis r 



It follows from ( pSf ) that ln(— ai) should be interpreted as ln(— ai + iO) = In(ai) + in. 
Besides, 

112(1 + ai) = 112(1 + ai - iO) = — ln(l + ai)[ln(ai) + iir] - Ll2(-ai) 

The result (^) is obtained for | A_l| ^ I Ay |. One can show that it remains valid in the case 
I Aj_| ~ |A||| if the expression (^ for A^ is used in (pS]). Actually, in eq. (|33| ) the difference 
between A and A^ is essential only in the overall factor 1/A^. For a screened Coulomb 
potential the amplitudes ( |5B| ) should be multiplied by the atomic form factor (1 — F(A^)). 
For the case of Moliere potential [|19| it reads 

1-^(A') = A^E^I^, (34) 

where 

«! = 0.1 , a2 = 0.55 , as = 0.35 , Pi = Pobi , (35) 

61 = 6 , 62 = 1-2 , h = 0.3 , po = mZ^l^jVll . 
Remind that the representation (PB|) is valid when |k2_L| , |ks_L| ^ m. 
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9 Cross section 



As it was suggested in [|T^], to overcome the problems of background in the measurement 
of photon sphtting one has to register the events with |f2,3| > /o where /o -C 1 is the angle 
determined by the experimental conditions. Let us consider the cross section integrated 
over fs for > /q. It is interesting to compare the exact result for this cross section 
{dcr/dx di2 ) following from (p3| ) and ( p6D with that obtained in the Weizsacker- Williams 
approximation {da approx / 2 )• The large logarithm corresponds to the contribution of 
the region A <^ p = |co'2f2 — ct;3f3)/2| where /a ~ xf2/{l — x). Taking the integral over A^ 
1% from A£„v 



m eq. 



up to Agjj, where 



see 



11) 



A' . = A 

mm 



icuj^x/2{l-x)f 



^eff 



and summing over the final photon polarizations we get for a pure Coulomb potential 
da, 



' approx 

dx dio 



SZ'^a^ g{x) /2(1 



/2 



X 



X 



-f2 — fo 



(36) 



For the case of a screened Coulomb potential the approximate cross section is 



da, 



approx 



dx dio 



AZ^a^ g{x) 



2 In 



'uJlXf2 

. (3o 



7 



X 



X 



-f2 — fo 



(37) 



The function 7 in eq. (0) is 

3 

7 = l-E«i(lna. + l)-2^ 



a,- In a,- 



aj In aj 



4 = 1 



i>j 



a,- — a, 



{31 



and the coefficients a^, 6j and /3o are defined in 



If A? 



JP^ > 1 then 7 



and eq. (|37| ) turns to eq. (|36D . If A 



'mm /Pi 

Coulomb potential the dependence of ctq ^da/dxdi2 on /2//0 is shown in Fig. 4 at /q = 10 



and X = 0.7 (curve 1), 



X 



2 < 1 then 7 = -0.15J 

fdxc 

0.3 (curve 2), where 



-ln(ALJ/3o^) 

For the case of a pure 

3 



and (^(x) is defined in (|3T|). For x = 0.7 the cross section given by (^) practically coincides 
with the curve 1. For x = 0.3 it is not the case (the curve 3 corresponds to the cross section 
in the Weizsacker- Williams approximation at x = 0.3). However, within a good accuracy 
the cross section da /dx at x = 0.3 agrees with that obtained from eq. (|36|). It should be 
so, since da/dx is invariant with respect to the substitution x — > 1 — x and at x = 0.7, as 
we pointed out above, the approximate result (|36|) is in accordance with the exact one. At 
X = 0.5 the difference between the exact result and the approximate one is rather essential 
(see Fig. 5). It can be explained as follows. The large logarithm appears as a result of 
integration with respect to fs over the range |(1 — x)f3 + xf2| -C x/2. After the integration 
over the azimuth angle ip between vectors £2 and —f^ we should integrate over from /q 
up to x/2/(l — x) and from x/2/ (1 — a;) to infinity. If x/2/(l — x) ~ fo then the contribution 
of the first region vanishes and the cross section becomes approximately two times smaller 
(in accordance with Fig. 5). 
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Since the amplitudes (|33D are obtained in zero-mass limit, it is interesting to estimate the 
accuracy of this approximation. Numerical calculations of the Born cross section da/dxdi2 
with the electron mass taken into account were performed in [jl^ (see Table V of that paper) 
for Z = 79, X = 0.87, uji = 1.7 GeV, 3.4 GeV, 6.1 GeV and for five values of the angle /2 
in the interval (1.2 — 2.8) mrad. These results are compared with ours in Table. One can 
see a good agreement everywhere. Only in one point corresponding to the smallest value 
of the transverse momentum k2± = 1.77 MeV the accuracy is 7%. 

For the case of a pure Coulomb potential the cross section da/dx can be approximated 
with a good accuracy by the following formula 



da, 



coul 



dx 



-f2 



^{x - 1/2) 



21n- 



2(1 



x] 



/o 



1 



+ (x ^ 1 



X 



(39) 



where 



Fix) 



1 X 

2 + — 



2x 



X 



-I ( 1 

xY \ X 



(40) 



If u;^/q//?q <C 1 then the corresponding expression for the cross section in a screened 
Coulomb potential reads 



da. 



dx 



^±^11A f2 In + 0.842 -F(x) 

x^ V /^o 



+ (x 1 — a;) 



(41) 



with F{x) defined in (^). The function F[x) characterizes the difference between the exact 
cross section and that obtained in the Weizsacker- Williams approximation. It is seen from 
Fig. 6 that this difference can be significant only at x being close to 0.5 . As for the total 
cross section, the difference between the exact result and the approximate one is about a 
few per cent. 

The inequality A ^ p which provides the applicability of the Weizsacker- Williams 
approximation corresponds to a small angle Lp between the vectors f2 and — fs ( when the 
vectors £2 and fs have almost opposite directions). So, it is interesting to consider the 
quantity da{ipmax) / dx which is the cross section integrated over the angle ip from —(pmax to 
ipmax- In the case of a pure Coulomb potential the dependence of {t^ fQ aQ)~^da{ipmax) / dx 
on ifmax is shown in Fig. 7 for different x and /o = 10"'^. One can see that the cross section 
becomes close to its total value at relatively large ^Pmax- The same conclusion is also valid 
for the case of a screened Coulomb potential. 
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Figure captions 



Fig. 1. Diagrams of the perturbation theory corresponding to the amphtude M^^\ eq. 

(!)■ 

Fig. 2. Diagrams of the perturbation theory corresponding to the amphtude M^'^\ eq. 

(I)- 

Fig. 3. Function g{x) from eq. (|30D for different polarizations: Qj, [x] (curve 1), 

(curve 2), (7+++(x) (curve 3), and g{x) (curve 4), eq. (pT]). 

Fig. 4. aQ^da/dxdf2 versus /2//0 for the case of a pure Coulomb potential, /o = 10^^, 
X = 0.7 (1), X = 0.3 (2), (To is given in the text. Curve 3 corresponds to the Weizsacker- 
Williams approximation for x = 0.3. 

Fig. 5. Same as Fig. 4 but for x = 0.5 (1). Curve 2 corresponds to the Weizsacker- 
Williams approximation. 

Fig. 6. Function F(x), eq. (|40|). 

Fig. 7. The dependence of {n (TQ)~^d(7{ipmax) / dx on (pmax for a pure Coulomb case at 
different x: x = 0.5 (1), x = 0.7 (2), and x = 0.9 (3); /o = lO'^. 
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Table. The photon-splitting cross section da/uidx in b/GeV sr for Z = 79, x = 0.87 



/2 

(mrad) 



cui = 1.7 GeV 
Present Paper [l7 
result 



ui = 3.4 GeV 
Present Paper ||T7 
result 



cji = 6.1 GeV 
Present Paper fLl 
result 



1.2 
1.6 
2.0 
2.4 
2.8 



22.72 
7.47 
3.09 
1.48 
0.793 



21.1 
7.4 
3.12 
1.51 
0.80 



3.11 
0.994 
0.404 
0.191 
0.101 



3.25 
1.03 
0.41 
0.19 
0.10 



0.565 
0.177 
0.0712 
0.0335 
0.0177 



0.56 
0.18 
0.072 
0.034 
0.018 
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